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ABSTRACT 



> 

iy-j | Steep soft X-ray (0.1-2 keV) quasars share several unusual properties: narrow 

£SJ , Balmer lines, strong Fell emission, large and fast X-ray variability, rather steep 2- 

■^j- ■ 10 keV spectrum. These intriguing objects have been suggested to be the analogs of 

CD ! Galactic black hole candidates in the high, soft state. We present here results from 

00 ! ASCA observations for two of these quasars: NAB0205+024 and PG1244+026. 

Both objects show similar variations (factor of ~ 2 in 10 ks), despite a factor of 
^ r-| about ten difference in the 0.5-10 keV luminosity (7.3 x 10 43 erg s _1 for PG1244+026 

O and 6.4 x 10 44 erg s _1 for NAB0205+024, assuming isotropic emission, Hq = 50.0 and 

qo = 0.0). 

The X-ray continuum of the two quasars flattens by 0.5-1 going from the 0.1-2 
keV band toward higher energies, strengthening recent results on another half dozen 



o 



I steep soft X-ray AGN. 

PG1244+026 shows a significant feature in the '1 keV region, which can be de- 
scribed by either as a broad emission line centered at 0.95 keV (quasar frame) or as 
k> \ edge or line absorption at 1.17 (1-22) keV. The line emission could be due to refiec- 

5_j ■ tion from an highly ionized accretion disk, in line with the view that steep soft X-ray 

quasars are emitting close to the Eddington luminosity. Photoelectric edge absorption 
or resonant line absorption could be produced by gas outflowing at a large velocity 
(0.3-0.6 c). 

Key words: Galaxies: Seyfert - Galaxies: individual: PG1244+026, NAB0205+024 
— X-rays: galaxies — Line: formation — Line: identification 



1 INTRODUCTION 

The ROSAT PSPC has found a large spread in the energy 
spectral indices of low-z quasars^]: 0.5 < ao.i-2teV < 3.5. In 
about 10% of cases Qo.i-2teV ~ 2 (e.g. Laor et al. 1994, 1997, 
Walter & Fink 1993, Fiore et al. 1994). The large spread in 
ceo.i-2keV favoured the discovery of its correlation with other 



properties. In fact, the steep soft X-ray quasars have then 
been realized to share a cluster of unusual properties: 

• narrow Balmer lines [j] ( Laor et al 1994, 1997, Boiler et 
al. 1995); 

• strong Fell emission (Laor et al 1994, 1997, Lawrence 
et al. 1997) 

• Rapid, large amplitude variability (factor of 2-50 on 



* We use "quasars" to describe broad line emission objects, re- 
gardless of luminosity. 



t the permitted lines have FWHM < 2000 km s _1 , yet still arc 
clearly broader than the forbidden lines. 
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timescales from minutes to months, Boiler et al., 1995, 
Brandt et al., 1995, Otani 1995, Boiler et al. 1997) 

• Somewhat steep hard X-ray spectra (2> ct2-iokeV > 
0.6, Pounds et al. 1995, Brandt et al., 1997); 

Pounds et al. (1995), suggest the latter to be a close 
physical analogy with the X-ray power-law produced by 
Comptonization in a hot accretion disk corona in Galac- 
tic black hole candidates (BHC) in their 'soft-high' state. 
This is not the only analogy between BHC and steep X- 
ray spectrum quasars. Laor et al. (1994, 1997) explained 
the correlation with H/3 FWHM as due to the larger size 
of a virialized broad emission line region for an AGN in a 
high L/LEdd state. Ebisawa (1991) found that while the soft 
component of 6 BHC observed by Ginga is roughly stable 
on time scales of 1 day or less, the hard component exhibits 
large variations down to msec time scales. These timescales 
translates to 10 4 years and 0.1 day for quasars, if they scale 
with the mass of the compact object. The soft component of 
BHC extends up to ~ 10 keV in BHC in 'soft-high' states, 
and it is often associated with optically thick emission from 
an accretion disk. If this is the case, the temperature should 
scale with the mass of the compact object as Mgj/ 4 , and 
the above energy translates to 0.1-0.4 keV for quasars. The 
rapid large amplitude variability shown by a few narrow line 
Seyfert 1 galaxies (NLSyl) at about 1 keV on timescales of 
hours to days (Otani 1995, Brandt et al. 1995, Boiler et al. 
1997) can then be analogous to the above BHC hard com- 
ponent flickering. 

A steep X-ray spectrum quasar with 10-100 times the 
luminosity of NLSyls, should be larger and so should vary no 
more rapidly than several days. Instead Fiore et al. (1998a) 
find that steep X-ray spectrum PG quasars commonly vary 
by a factor 2 in 1 day. Variability seems therefore correlated 
with X-ray spectral slope and Balmer line width (and there- 
fore possibly with the accretion rate) rather than with the 
luminosity. 

Evidence for spectral features in the '1 keV region in 
many steep soft X-ray quasars is building up (Turner et al., 
1991, Brandt et al., 1994, Otani et al., 1995, Comastri et al., 
1995, Leighly et al., 1997, 1998a,b). Instead, 'normal' Seyfert 

1 galaxies (having broad Balmer lines and flatter soft X- 
ray spectra) usually have their strongest absorption features 
at lower energies (in the 0.6-0.9 keV 'oxygen' band). An 
intriguing possibility is that the appearance of these features 
at different energies also depend on L/LEdd- 

Detailed high energy X-ray spectra of luminous quasars 
with steep soft X-ray spectra are essential to understand 
the 'narrow-broad line' phenomenon in AGN, in particular 
whether the peculiar X-ray properties depend on optical lu- 
minosity, optical-to- X-ray ratio (aox), or on their Edding- 
ton ratio. To this end we selected two bright quasars with 
ao.i-2fav >2.0 (Fiore et al., 1994) at the extreme values of 
optical luminosity, both with low Galactic Nh (Table 1) of 
1.9 x 10 20 cm" 2 for PG 1244+026, and of 3.0 x 10 20 cm" 2 for 
NAB0205+024, Elvis et al., 1989) and observed them with 
ASCA. We report the results in this paper. 

2 OBSERVATION AND DATA REDUCTION 

Table 1 gives the redshift, My, the 0.2-2 keV luminosity, 
aox, the average PSPC count rate and spectral index and 



the Galactic Nh for the two quasars. Table 2 gives the ASCA 
observation log, the SIS and GIS exposure times and count 
rates. 

Both observations were performed in two CCD mode 
with the source at the '1CCD mode' position. Data reduc- 
tion was performed using FTOOLS 3.6. We used "bright" 
mode SIS data, combining LOW, MEDIUM and HIGH bit 
rates. Conservative cleaning criteria were applied (minimum 
Earth occultation = 7 degrees, minimum magnetic rigid- 
ity = 6 GeV/c, minimum bright Earth angle = 20 degrees, 
and excluding data collected in the first 32 seconds after 
the satellite passage in the SAA and through the day-night 
terminator). Counts, light curves and spectra from the two 
quasars were accumulated in circular regions of 3 and 4 ar- 
cmin radius for SIS and GIS respectively. 

We are interested in the high energy spectrum of these 
sources and since they are rather faint, and possibly very 
steep, background subtraction plays a crucial role. Back- 
ground counts were accumulated from regions surrounding 
the sources and compared with counts accumulated from the 
same regions from 'blanksky' observations. The 'local' and 
'blanksky' background counts were always within 10 % of 
each other for the four ASCA instruments. To obtain the 
best possible signal to noise in the background subtracted 
spectra we therefore used the 'blanksky' background in our 
spectral analysis. We extracted background spectra from 
'blanksky' event files using the same regions as for source 
extraction. The count rates of the two sources become that 
of the background at about 7 keV (observer frame). After 
background subtraction PG1244+026 is observed in the GIS 
up to 10 keV and NAB0205+024 up to 8 keV (9.3 keV quasar 
frame), both at the > 3 a level. 

Spectral fits were made separately to the spectra from 
the four ASCA instruments and to the spectra obtained 
combining together the data from the two SIS and GIS de- 
tectors. The results were consistent with each other. In the 
following we present the results obtained following the sec- 
ond approach. In some cases x 2 are smaller than 1. This 
is due to the prescription adopted in adding the spectra, 
for the propagation of the errors (the Gehrels, 1986, algo- 
rithm: error = 1.0 + SQRT(N + 0.75)). Spectra were always 
rebinned following 2 criteria: a) to sample the energy reso- 
lution of the detectors with four channels at all energies 
where possible, and b) to obtain at least 20 counts per en- 
ergy channel. In the spectral fits we limited ourselves to the 
0.6-10 keV energy band, to minimize the systematic effect 
due to the uncertainty in the SIS calibration below 0.6 keV. 
Nh was always constrained to be greater than or equal to 
the Galactic value along the line of sight. 

In all cases the quoted errors represent the 90 % confi- 
dence interval for 1 interesting parameter. 



3 VARIABILITY 

ASCA observed the two quasars for a total elapsed time of 
105 ks (PG1244+026) and 119 ks (NAB0205+024). We can 
therefore study the variability of these sources on time scales 
from a few hundred seconds to about 1.5 days. 

Figure 1 and 2 show the SIS light curves of the two 
sources. Following Ptak et al. (1994) the values were com- 
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Table 1. Steep soft X-ray quasars 



name 


z 


M« 


Lo.2-2keV 




PSPC 


a 0.1-2feeV 










10 45 erg s" 1 




cts/s 




10 20 cm- 2 


PG1244+026 


0.048 


-21.1 


0.14 


1.4 


1.0 


2.3 ±0.1 


1.9 6 


NAB0205+024 


0.155 


-25.0 


1.8 


1.6 


0.7 


2.3 ±0.1 


3.0 6 



a H = 50, <?o = 0; 6 Elvis et al. 1989 



Table 2. ASCA observations 



name 


Dates 


Exposure SIS-GIS 


2 SIS count rate 


2 GIS count rate 






ks 


cts/s 


cts/s 


PG1244+026 


1-3 Jul 1996 


37-39 


0.442+/-0.004 


0.241+/-0.004 


NAB0205+024 


18-20 Jan 1996 


50-54 


0.150+/-0.001 


0.153+/-0.002 



D PG1244+026 
d I I 1 




w o 



2x10* 4xl0 4 6xl0 4 8x10* 10 5 

Elapsed Time (s) 

Figure 1. The SIS 0.5-7 keV light curve of PG1244+026. 



NAB0205+024 




5x10* 10 s 

Elapsed Time (s) 

Figure 2. The SIS 0.5-7 keV light curve of NAB0205+024 

puted using variable bin sizes corresponding to good time 
intervals longer than 200 seconds. 

Significant variations are evident in both light curves. A 
factor of two variability, both down and up, in about 10 ks is 
present at the beginning of the light curve of PG1244+026. 



A drop of a factor of about 2 is present at the beginning of 
the light curve of NAB0205±024 on a ~ 15 ks timescale. In 
this paper we limit ourselves to pointing out that roughly 
similar variability is observed in two sources which differs 
in luminosity by a factor ~ 10. Significant, although rather 
small, spectral variability is also present in the ASCA ob- 
servations of the two sources. A systematic analysis of the 
variability in different energy bands and of the spectral vari- 
ability is in progress and will be presented, together with a 
similar analysis on a sample of about 20 Seyfert 1 galaxies 
and quasars observed by ASCA and BeppoSAX, in a paper 
in preparation. In the following sections we present the aver- 
age properties of the spectra. We anticipate that the spectral 
variability will not modify the results presented here. 



4 SPECTRAL ANALYSIS 

Figures 3 and 4 show the SIS±GIS spectra of PG1244+026 
and NAB0205+024 fitted with a simple power law absorbed 
at low energy by a column of cold gas equal to or higher 
than the Galactic column along the line of sight. Table 3 
gives the best fit parameters and the x 2 . It is clear that 
this simple model is inadequate to describe the 0.6-10 keV 
spectrum of both quasars. A hard tail, larger than the 5-10 
% systematic uncertainties at these energies (Gendreau & 
Yaqoob 1997), is evident in both cases. In the spectrum of 
PG1244+026 there is also a significant excess with respect 
to the model about 1 keV. The feature is visible in both SIS 
and GIS detectors. We discuss these two findings in turn in 
the next two sections. 



4.1 0.6-10 keV continuum 

The ASCA SIS and GIS responses are peaked at 1-2 keV 
and decrease sharply at higher energies. Therefore the ASCA 
'2-10 keV slopes are strongly biased toward the lowest en- 
ergy boundary. This means that some caution should be 
used when comparing ASCA '2-10 keV slopes with those 
of experiments whose responses peak around 6 keV such 
as EXOSAT, GINGA, BeppoSAX and XTE. To address 
this, we fitted the SIS±GIS spectra of PG1244+026 and 
NAB0205+024 with a simple power law in the observed 0.6- 
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Table 3. PG1244+026 & NAB0205+024 spectral fits 



model 



N? 



a E or T b a H or tt/2n 



ore 



or A c 



x 2 (dof) 



PG1244+026 

PL 0.6-10 keV 

PL 2-10 keV 

PL 3-10 keV 

PL 4-10 keV 

Broken PL 0.6-10 keV 

PL+Raym 0.6-10 keV 

PL+BB 0.6-10 keV 

PL+Comp.Refl 0.6-10 keV 



1.9+0.6 
1.9FIXED 
1.9FIXED 
1.9FIXED 
3.6l 2 ; 6 7 
1.9+0.7 

c o+8.0 
°-°-3.4 
2 6+ 1 - 3 



1.67+0.04 
1.35+0.12 
1.03+0.30 
0.67+0.55 
1.80lJ;i 
0-88l°- 2 5 
0.16+0.03 
2.78+0.06 



1.06"* 



„+0.24 
-0.38 

1.54+0.06 
1.40+0.15 

> 8 



2.8+0.6 
> 0.5 



177.9 (148) 
48.4 (82) 
29.1 (58) 
14.4 (38) 

149.1 (145) 
135.3 (145) 

139.2 (145) 
150.5 (146) 



NAB0205+024 

PL 0.6-10 keV 3.0+0.1 

PL 2-10 keV 3.0FIXED 

PL 3-10 keV 3.0FIXED 

PL 4-10 keV 3.0FIXED 

Broken PL 0.6-10 keV 3.0+1.4 

PL+Raym 0.6-10 keV 3.0+7.0 

PL+BB 0.6-10 keV 3.0+0.2 
PL+Comp.Refl 0.6-10 keV 3.0+0.2 



1.38+0.02 
1.09+0.10 
0.94+0.20 
0.62+0.45 
1.5012;? 
0.49+0.05 
0.16+0.02 
1.44+0.04 



0.98+0.10 
1.02+0.25 
1.15+0.07 
> 4.2 



2.4+0.5 
< 0.01 



126.2 (147) 
42.4 (86) 
23.3 (58) 
9.5 (38) 
85.2 (145) 
82.8 (144) 
86.2 (145) 
95.2 (146) 



10 20 



in keV; c metal abundances 



PG1 244+026 



NAB0205+024 



ASCA SIS - CIS 




ASCA SIS - CIS 




channel energy (keV) 



channel energy (keV) 



Figure 3. The SIS+GIS spectra of PG1244+026 fitted with a 
simple absorbed power law model. The lower panel shows the 
ratio between the data and the best fit model. 



Figure 4. The SIS+GIS spectra of NAB0205+024 fitted with 
a simple absorbed power law model. The lower panel shows the 
ratio between the data and the best fit model. 



10, 2-10, 3-10 and 4-10 keV ranges. The results are given in 
Table 3 and show in Figure 5. 

In the last three series of fits the Nh was fixed to the 
Galactic value. The best fit slopes flatten by 0.5-1 going 
from the low energy dominated to the higher energy range. 
(We note that for NAB0205+024 the quasar redshift implies 
that these slopes refer to slightly harder energy ranges: 10 
keV in the observer frame corresponds to 11.6 keV in the 
quasar frame). 

The PSPC spectral index for both quasars is 2.3(±0.1, 
Fiore et al., 1994), steeper than any of the values in table 
3. For NAB0205+024 the ASCA low energy index is flatter 
than the PSPC by 0.4-0.8, and the ASCA best fit N H coin- 



cides with the Galactic value, suggesting that the spectrum 
continues to steepen below ASCA X-ray energies. 

To parameterize a curved spectrum we have also fitted 
the 0.6-10 keV spectra of the two quasars with: a broken 
power law model; a power law + optically thin plasma emis- 
sion (Raymond & Smith 1977) model; a power law + black 
body model; a power law + Compton reflection model. The 
results are again in Table 3 and confirm the presence of sig- 
nificant curvature in the spectra of these quasars. 

In PG1244+026 the power law + Raymond-Smith 
model gives a \ 2 significantly better than the broken power 
law model because it fits also the ~ 1 keV feature (see next 
section). However, the best fit energy index is still very steep 
and positive residuals are evident above 4 keV (Figure 6). 
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Figure 5. The PG1244+026 (filled circles) and NAB0205+024 
(open triangles) best fit power law index in the observed 0.6-10, 2- 
10, 3-10 and 4-10 kcV ranges as a function of the range minimum 
energy. The NAB0205+024 points have been shifted from real 
energy for a sake of clarity. The PSPC 0.1-2 keV points are also 
showed for comparison. 



In NAB0205+024 the power law + Raymond-Smith model 
gives an acceptable \ 2 only for low metal abundances, given 
the absence of significant features in the spectrum, similar 
to the ROSAT results of Fiore et al. (1994). Above 4 keV the 
best fit is similar to that in the broken power law model. Fits 
with a power law + black body models give acceptable \ 2 
in both cases. The power law indices in this case are slightly 
steeper than in the broken power law model. Inspection of 
the residuals shows again a slight excess of counts at high 
energy. In PG1244+026 the power law + black body model 
again gives a \ 2 significantly better than the broken power 
law model because it partly fits also the feature around ~ 1 
keV (see next section) . Fits with a Compton reflection model 
(plrefl in xspec) give x 2 significantly higher than the pre- 
vious models and push the parameter Q/2-k to high and 
implausible values. 

We do not see any significant line emission at the en- 
ergies of the iron Ka lines. The 90 % upper limits to the 
equivalent width of a narrow line at 6.4 (and 6.7 keV), rest 
frame, in PG1244+026 and NAB0205+024 are 400 eV and 
230 eV (640 eV, 314 eV) respectively. 

The x 2 for the fits to the PG1244+026 spectra are much 
higher than those for the similar fits to the NAB0205+024 
spectra because of the presence of the T keV feature in the 
former quasar. We discuss this feature next. 



PG1 244+026 

ASCA SIS - CIS 




1 2 
channel energy (keV) 



Figure 7. The SIS+GIS spectra of PG1244+026 fitted in the 
0.6-4 keV range. The lower panel shows the residuals when the 
spectrum is fitted with a simple absorbed power law 



4.2 Low energy features 

We have already noted that a '1 keV feature seems to be 
present only in PG1244+026, the quasar with the lower lu- 
minosity. Figure 7 shows the results of a simple power law 
model fit to the 0.6-4 keV spectrum of PG1244+026: the '1 
keV feature is clearly visible. 

It has been recently realized that Residual Dark Distri- 
bution (RDD, the error in ASCA's onboard correction for 
CCD dark current) can affect the low energy ASCA SIS 
spectra. The symptom of the RDD problem is a sudden de- 
crease in the SIS effective area towards lower energies from 
about 1 keV. We are however confident that this has little 
effect on the '1 keV feature for the following reasons: a) SIS 
data of PG1244+026 and NAB0205+024 have similar RDD 
value but the '1 keV feature is visible in the spectrum of 



© 0000 RAS, MNRAS 000, 000-000 



6 F. Fiore et al. 



PG1244+026 only, the 90 % upper limit to any unresolved 
gaussian line emission at about 1 keV in NAB0205+024 is 
of 20 eV; b) the feature is present also in the GIS spectrum; 
the feature is present at the same level in SISO and SIS1, 
while the RDD tends to be greater for SIS1 than for SISO. 

To investigate the physical origin of the '1 keV feature 
and to quantify its strength we performed a series of spectral 
fits to the SIS and GIS spectra of PG1244+026 using only 
the 0.6-4 keV range, to avoid the complications of the hard 
tail (see previous section). The results are reported in Table 
4, along with those of the other fits in this section. 

4-2.1 Emission Line model 

The fit with a power law model plus a gaussian line gives a 
small x 2 (90.1 for 101 degrees of freedom). The equivalent 
width of the 1 keV feature, in the power law plus gaussian 
line fit, is similar in the SIS and GIS detectors. In Table 4 we 
report the SIS determination of 64±15 eV. The line width 
is well constrained and cleanly resolved in the SIS spectrum 
to a = 0.11±g;g| keV. 

A fit with a power law plus a thermal plasma model 
(Raymond-Smith 1977) gives an acceptable x 2 (93.1, 100 
dof). The best fit temperature (1 keV) implies that the '1 
keV emission is dominated by a blend of iron-L and neon 
emission lines. 

An emission line feature can be mimicked by fitting a 
spectrum with a strong absorption feature at slightly higher 
energies. We have then fitted the SIS and GIS spectra with 
models including absorption structures. 



4.2.2 Absorption Edge Fits 

Fits with one absorption edges gives \ significantly higher 
than the previous case (102.5, 102 dof). The edge energy 
(1.17 keV, quasar frame) is consistent with that of NelX 
and/or iron L FeXVI and FeXVII. The best fit neutral N H 
is significantly higher than the Galactic value. This is reason- 
able, since if there is highly ionized Ne and Fe L absorption 
it is likely to have also highly ionized oxygen absorption at 
0.74-0.87 keV. We then refitted the SIS and GIS spectra with 
a model including three absorption edges, at the energies of 
the most abundant ions in highly ionized gas with high NelX 
abundance: OVII, OVIII and NelX-FeXVI-FeXVII (Nicas- 
tro et al. 1998), fixing the cold Nh to the galactic value. The 
results were not satisfactory. The depth of the oxygen edges 
is zero with small upper limits and the \ 2 ls significantly 
higher than in the previous case: 114.7. Leaving Nh free 
improves the x 2 \ but the oxygen edge depths are still zero 
and the fit resembles completely the single edge fit. Fixing 
Nh to the Galactic value but leaving free the energies of 
two edges produces again a good fit (x2 = 94.8, 101 dof). 
The best fit energy of one edge is again 1.12 keV, but that 
of the other edge is < 0.64 keV (observer frame), close to 
the lower boundary of the observed range. So, if the cold 
absorption is fixed to the Galactic value, then there must 
be additional absorption edge(s) at energies lower than the 
observed range, corresponding to oxygen less ionized than 
OVI. We note however that this conclusion is weakened by 
the unknown contribution of the SIS RDD, which pushes 
low energy events below the detection threshold. 



A 1.17 keV absorption feature can also be interpreted in 
terms of blueshifted oxygen absorption (Leighly et al. 1997). 
In this case, assuming that the absorption is mostly due to 
OVIII, the shift from the quasar frame would be equivalent 
to z=-0.38. A more complex continuum has little effects on 
the best fit parameters of absorption edges. 



4-2.3 Ionized absorber models Fits 

We fitted the data with a detailed ionized absorber model 
(not including resonant scattering absorption lines) . We first 
generated a grid of photoionization equilibrium models us- 
ing Cloudy (Fcrland 1996), and fitted the spectrum in- 
terpolating by this grid, using the method of Fiore et al. 
(1993). To calculate the models we have assumed the ob- 
served spectral energy distribution (Fiore et al. 1995, Elvis 
et al. 1994). This is important, since the soft X-ray spectrum 
of this source strongly differs from that of 'normal' Seyfert 1 
galaxies, where warm absorbers are usually found (Reynolds 
1997). A steep soft X-ray spectrum can completely ionize 
oxygen and neon but not iron, and so can produce edges in 
the 1-2 keV (Fe-L) and 7-9 keV (Fe-K) ranges, but not in the 
'oxygen' 0.6-0.9 keV band. A fit with this model produces 
an acceptable \ 2 ( see Table 4, ionized absorber model 1). 

In Figure 8 we show the best fit steep SED model (thick 
line), and a photoionization model obtained using a stan- 
dard, much flatter AGN SED (a power law of a = 1.2 from 
UV to X-rays, thin line), which, above 1 keV, gives a compa- 
rably good fit to the data. While iron in the flat SED model 
has a ionization structure similar to that of the model ob- 
tained using the right SED, oxygen is much less ionized: note 
the deep OVIII edge present in the flat SED model. 

In this fit the redshift of the absorber is significantly 
higher than that of the quasar, because the main feature in 
the transmitted spectrum is the FeXVIII edge at 1.36 keV, 
while the deepest edge in the quasar spectrum is at 1.17 keV 
(quasar frame). However, a good fit can be also obtained for 
a different absorber redshift (z=-0.33, Table 4, ionized ab- 
sorber model 2). In this case the 1.17 feature is interpreted 
in terms of OVII and OVIII absorption. We cannot discrim- 
inate between these two solutions on statistical grounds. 

We have also tried fits with a collisional equilibrium 
model (Nicastro et al. 1998). The results were very similar 
to those obtained in the case of photoionization equilibrium 
(see Table 4, ionized absorber model 3). 

In all fits with detailed warm absorber models the col- 
umn of cold gas is significantly higher than the Galactic, 
value, similar to the values found using single edges to pa- 
rameterize the absorber (Table 4). 



4.2.4 Absorption Line Fits 

Fits with a single gaussian absorption line do not give ac- 
ceptable \ 2 - Fits with two or more absorption gaussian lines 
can produce \ 2 of 96 or smaller. These models are indistin- 
guishable for the SIS from models with a broad absorption 
notch, which we discuss in the following. Physical models 
including resonant absorption lines, as well as absorption 
edges from ionized plasma will be discussed in a paper in 
preparation (Nicastro et al. 1998b). 



© 0000 RAS, MNRAS 000, 000-000 



ASCA observations of two steep soft X-ray quasars 7 



Table 4. PG1244+026: spectral fits in the band 0.6-4 keV 



Continuum models 
model 


N H 


rpb 


A c 




X ( aot ) 


PL 

PL + BB 
PL + Raym 


1.70±0.02 
1.50+0.11 
1.61+0.03 


1.9+1.0 

Q l^L — 1— n 

y .ozco.u 
1.9+0.3 


n i a -4- n no 
1.00+0.06 


>0.5 


- 
- 


131.8 (104) 

1 no ^ (~\ ni ~i 
±uy.o ^ iui ) 

93.1 (100) 


emission line 






line K b 


11I1U J_j V V 




X ( ao1 ) 


PL + Gauss 


1.60+0.04 


1.9+0.5 


0.91+0.03 


64+15 


- 


90.1 (101) 


absorption edge 


models 




edge E b 


T 




X 2 (dof) 


PL + 1 edge 
PL + 2 edges 


1.88+0.11 
1.88+0.11 


8.2+3.0 
1.9FIXED 


1.12+0.03 
1.12+0.10 
< 0.64 


0.25+0.07 
0.26+0.08 
0.51+0.14 




102.5 (102) 
94.8 (101) 


ionized absorber models 




logA^ (ionized) 


U or logT 


/ abs. 


X 2 (dof) 


1 - Phot. ion. 

2 - Phot. ion. 

3 - Coll. ion. 


1.92+0.10 
1.94+0.10 
1.90+0.10 


11.1+3.5 
11.0+3.5 
11.3+3.5 


22.59+0.15 
21.56+0.10 
22.62+0.14 


2.00+0.12 
1.20+0.10 
6.92+0.07 


0.25+0.07 
-0.33+0.10 
0.24+0.06 


93.2 (101) 
95.7 (101) 
90.7 (101) 


absorption notch models 




T b , notch E b 


notch width d 


Cov. frac. 


X 2 (dof) 


PL + notch 
PL+BB+notch 


1.93+0.08 
1.48±0.20 


9.3+3.3 
7.8+4.4 


1.44+0.05 
T = 0.16 ±0.03 
1.16+0.03 


680+70 
14+7 


0.14+0.04 
1FIXED 


97.5 (100) 
97.1 (100) 



a in 10 20 cm 2 ; 6 in keV; c metal abundances; d in eV 



4-2.5 Absorption Notch Fits 

Fits with an absorption notch give \ 2 higher than those 
with an emission lines by A\ 2 w 7. While the power law + 
notch fit is formally acceptable, the best fit value of the notch 
width is implausibly large (almost 1 keV), forced by the very 
low value of the covering fraction required by the fit. How- 
ever, the notch best fit parameters are strongly dependent on 
the proper modeling of the continuum. For example, using 
a power law + a black body to parameterize the continuum 
gives an acceptable fit fixing the notch covering fraction to 
1, which in turn results in a much more reasonable value for 
the notch width of 14+7 eV. The ASCA band width in not 
wide enough and its spectral resolution is not good enough 
to constrain adequately both a complex continuum and the 
notch parameters. 



4-2.6 Comparison with the PSPC results 

PG1244+026 was observed with the PSPC in December 
1991 and the results of this observation have been reported 
by Fiorc ct al. (1994). The 0.6-2 keV flux level during the 
PSPC observation was about 50 % lower than the mean flux 
in the ASCA observation. The fit of an absorbed power law 
to the PSPC spectrum gives an acceptable \ 2 (22.8 for 25 
dof), a .i-2fcef = 2.3±0.1 and N H = 2.9± 0.3 x 10 20 cm" 2 , 
slightly higher than the Galactic value. Most of the PSPC 
counts were detected below 0.3 keV, in the 'Carbon' band. 
The quality of the spectrum between 0.3 and 2 keV is not 
very high and emission or absorption features fainter than 



~ 20% cannot be excluded in this energy band. No evidence 
of spectral variability is present in the ROSAT data despite 
a factor of 2 flux variability. 

The PSPC data strongly constrain the level of any cold 
or warm absorption affecting the 'Carbon' band. Best fitting 
models to the ASCA data including absorption features in 
the 1-2 keV band (Table 4) require a rather large absorp- 
tion in addition to the Galactic one below 1 keV. Therefore, 
it is important to study whether the ASCA best fit models 
are consistent with the PSPC ones. Rather than performing 
joint fits to the ASCA and PSPC data, which are compli- 
cated by the large uncertainty in the relative PSPC/ASCA 
SIS calibration, and by the detailed shape of the continuum 
over the broad 0.1-4 keV band, we fitted the PSPC data 
with a power law model including the emission or absorp- 
tion features found in the previous section (see Table 4). The 
results are in Table 5, where we also report (in brackets) the 
99 % parameter upper limits, or confidence intervals, when 
appropriate. 

We see that the presence of an emission line at 0.91 keV 
is not required by the PSPC spectrum, but an equivalent 
width of 65 eV is not excluded (10 % probability). However, 
the presence of a cold absorber of thickness 8.2 ± 3.0 x 10 2t) 
cm" 2 is inconsistent with the PSPC spectrum (probability 
< 1%), while the presence of an edge at 0.62 keV with r = 
0.51 ± 0.14 is only marginally consistent with the PSPC 
result. 
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Table 5. PG1244+026: spectral fits to the PSPC 



model 


N a 


OLE 


line or edge E 6 


line EW C or r 


X 2 (dof) 


PL 

PL + gauss 
PL + 1 edge 
PL + 2 edges 


2.9±0.3 
2.9±0.3 
2.7±0.4(0.6) 
2.7±0.4 


2.3±0.1 
2.3±0.1 
2.2±0.2 
2.2±0.2 


0.91FIXED 
1.12FIXED 
1.12FIXED 
0.62FIXED 


0+68 (+100) 
3+ ' 4 
3+ - 4 
0+0.34(+0.62) 


22.8 (25) 

22.8 (23) 

20.9 (24) 
20.9 (23) 




a in 


10 20 cm- 


2 ; 6 in kcV; c in cV 








Energy (keV) 



Figure 8. The thick line shows the best fit photoionized absorber 
model for PG1244+026 using the observed steep spectrum (note 
that the absorber has z=0.25, while the quasar has z=0.048). The 
thin line shows a photoionization model obtained using a much 
flatter SED (see text). Both models produce Fc-L absorption fea- 
tures around 1 keV. The flat SED also produces deep OVIII edge, 
which is not observed 

5 DISCUSSION 
5.1 Continuum 

The ASCA observations of PG1244+026 and NAB0205+024 
have shown that the X-ray continuum of these two quasars 
flattens by 0.5-1 passing from the 0.1-2 keV (PSPC) to the 
2-10 keV band. Similar results were obtained by Brandt, 
Mathur & Elvis (1997); and by Comastri et al. (1998) and 
Leighly et al. (1998a) on TONS180, Pounds et al. (1995) and 
Fiore et al (1998b) on REJ1034+390, Leighly et al. (1998b) 
on AKN564. It appears that the X-ray spectrum of a sizeable 
number of steep PSPC and narrow Balmer line quasars has 
significant curvature, being flatter at higher energies. 

This could be due to different components influencing 
the spectrum at different energies, as might happen in 'nor- 
mal' broad lines Seyfert 1 galaxies and quasars, where a soft 
component is often present. The relative intensity of the two 



components would be quite different from 'normal' quasars. 
Laor et al. (1997) suggested a fainter hard component rela- 
tive to the optical in the majority of low redshift PC quasars 
(assuming a two component model). However, Grupe (1996) 
found evidence for a stronger soft excess in a sample of soft 
X-ray selected Seyferts dominated by narrow-line objects. 

A large relative intensity of the soft component has im- 
portant consequences on various competing models for the 
soft component. Disc reprocessing models (Matt et al. 1993, 
Fiore et al. 1997) would require highly anisotropic emission 
to account for the discrepancy between the observed soft 
and hard fluxes. Optically thin free-free emission (e.g. Bar- 
vainis 1993) is ruled out by these observations because the 
best fit power law slope is still too steep to fit the spec- 
trum above 4 keV, because variability rules out optically 
thin plasma (Elvis et al. 1991) and because of the implausi- 
bly low metal abundances (< 1% solar, see Table 3) required 
in NAB0205+024 (see Sect. 3.2.1). The most likely origin for 
the steep component is Comptonized disc emission (e.g. Cz- 
erny & Elvis 1987, Fiore et al. 1995, Pounds et al. 1995). 

The high energy spectral index of the two quasars an ~ 
1.0 (see Table 3) is consistent with that of 'normal' Seyfert 
1 galaxies (e.g. Nandra & Pounds 1994, Nandra et al. 1997). 
The error on o>h is however large and so no strong conclusion 
can be drawn on the origin of the hard emission. An answer 
to this question must await the large area and high energy 
sensitivity of AXAF, XMM and Spectrum X-gamma. 

5.2 Origin of the ~ 1 keV feature in PG1244+026 

The ~ 1 keV feature in PG1244+026 could be explained 
in terms of either a broad (a = 0.1 keV) emission line at 
0.91 keV (0.95 keV quasar frame) of about 60 eV equivalent 
width or a r = 0.25 absorption edge at 1.17 keV (or an 
absorption notch at 1.22 keV). These possibilities cannot be 
discriminated between on statistical grounds. 

An absorption interpretation requires additional low en- 
ergy absorption, either cold (with a column density higher 
than Galactic by ~ 7 x 10 20 cm~ 2 (see Tables 4 and 5), in 
contrast with the PSPC results, or a warm absorber with a 
peculiar ion abundance distribution. 

The observed spectrum can be interpreted in terms of 
either an inflowing (v/c=0.25) or an outflowing (v/c=-0.33) 
absorber. In the first case the ion contributing most to the 
absorption is FeXVIII, in the second case it is OVIII. The 
two cases cannot be discriminated on statistical grounds. 
While an outflowing absorber has been suggested in several 
other cases (see e.g. Mathur et al. 1994), this would be the 
first case for an inflowing highly ionized absorber. 

A similar situation is found in IRAS 13224-3809 by 
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Otani et al. (1995), in AKN564 by Brandt et al. (1994) and 
by Leighly et al. (1998a). Otani et al (1995) and Leighly et 
al. (1997) interpret the features in the 1-2 keV band in terms 
of blueshifted absorption from relativistically (v=0.2-0.6 c) 
outflowing material. If the 1.17 keV feature in PG1244+026 
is due to a blueshifted OVIII photoelectric absorption, then 
the absorption seen below 0.64 keV may be due to CVI pho- 
toelectric absorption from the same gas. This distribution is 
far from an equilibrium distribution (see eg. Nicastro et al. 
1998), not an impossible situation considering the large vari- 
ability observed in this source. 

The absorption features seen between 1 and 2 keV may 
also be interpreted in terms of resonant lines (e.g. Leighly 
et al. 1997). If the ion producing the absorption is oxygen 
OVIII (resonant line at 0.65 keV), the best fit notch energy 
of 1.22 keV (quasar frame) implies a very high gas velocity: 
-0.56c. Fe XVIII (E=0.87 keV) or Fe XVII (E=0.81 keV) 
can also contribute to the absorption, because of their high 
oscillator strengths, 1.7, 0.6 respectively (Kato et al. 1976), 
and abundances. If the 1.22 keV absorption notch is due to 
these ions then the velocity of the outflowing gas will be 
smaller, -0.33c. 

In any case, in the blueshifted absorption scenario the 
gas is outflowing at velocities which are a sizeable fraction 
of c, reminiscent of blobs of gas in jets. It is interesting to 
note that similar absorption features have been observed in 
Blazars (e.g. PKS2155-304, Canizares & Kruper, 1884, other 
BILacs, Madejski et al. 1991, 3C273, Grandi et al. 1997), but 
usually below 1 keV. Somewhat surprisingly this implies less 
extreme conditions in these radio-loud objects than in our 
radio-quiet quasars. High redshift radio loud quasars may 
have similar jet-related absorption too (Elvis et al., 1997). 

An alternative interpretation of the 1 keV feature is in 
terms of an emission line due to highly ionized oxygen, neon 
(NelX) and/or to Iron L. There are two possible origin for 
this line: recombination in an optically thin thermal plasma, 
or "reflection" in photoionized matter. 

5.2.1 Thermal plasma 

A thermal plasma is highly implausible on physical grounds 
(Elvis et al. 1991, Fiore et al. 1995). Emission measure is 
~ 2.5xl0 65 cm~ 3 , For a spherical source (with radius R) 
and constant electron density n e , n e = 2.43 x 10 32 -R~2 
cm -3 . R < 10 14 cm from the observed X-ray variability, 
n e < 1.5xl0 12 cm~ 3 , This implies an electron scattering 
optical depth tt > 30. With such values a thermal plasma 
is no longer optically thin. The situation is even worse if 
the matter is clumpy, as the density of each cloud must be 
greater. 

5.2.2 Photoionized matter 

The second possibility is that the '1 keV emitting matter 
is photoionized by the central nucleus. We can assume that 
the gas is not covering the source because there is not signif- 
icant OVII and/or OVIII absorption. We therefore assume 
that we are not observing a "warm absorber" , but rather 
a "warm reflector". This could either be a warm absorber 
viewed from its side, in which case the matter would be opti- 
cally thin to Thomson scattering; or the accretion disc, and 



the matter would be thick. Since both the "reflector" and 
the primary emission (which provides most of the contin- 
uum) are observed, the optically thick case gives the highest 
values of the equivalent width. However, even in the opti- 
cally thick case the expected EWs can barely account for the 
observed values (see below); hence we neglect the optically 
thin case altogether. 

The observed line (which is significantly broad) may be 
a blend of several lines (see for instance Zycki et al. 1994, 
Netzer 1997). The observed energy suggests the 0.92 keV Ne 
IX and 1.02 keV Ne x recombination lines, the 0.87 keV Oxy- 
gen viii recombination (to ground state) line and the iron 
L (around 0.8 keV) lines being the most important. None 
of these lines alone can account for the observed EW: for 
instance, the maximum value (i.e. for a face-on disk with a 
intervening ion fraction ~ 0.6) for the Ne IX line is about 10 
cV, while that of the oxygen recombination line is about 15 
eV (note that a O viii Ka recombination line at 0.65 keV 
with a similar EW should also be present; the 90% upper 
limit on such a line is 30 eV). These values have been cal- 
culated using the formulae of Basko (1978), and assuming a 
reasonable ionization structure. (In the disc hypothesis an 
iron Ka line at 6.5-6.9 keV is also expected, but the up- 
per limit of 300-400 eV does not exclude the presence of 
such a line.) Allowing for a possible factor of 2 neon and/or 
iron overabundance (an oxygen overabundance would de- 
crease the Ne line while not increasing the O line) and/or 
anisotropy of the illuminating radiation, the observed equiv- 
alent width could be explained (note that in many Seyfert 
1 galaxies the iron Ka line is also stronger than expected, 
suggesting iron overabundance or anisotropic illumination) . 

In this scenario the '1 keV feature may arise from an 
highly ionized accretion disk. In the Matt et al. (1993) mod- 
els high ionization is mainly due to a high accretion rates 
(the ionization parameter depends on rh 3 ). The detection 
of these emission lines in steep soft X-ray quasars would 
then be further evidence of high L/Z/Edd- We note that since 
recombination can occur only in highly ionized atoms, we 
would not expect features of this kind in 'normal' quasars, 
as their disc should be much less ionized, as in fact observed. 

The strong dependence of the ionization parameter on 
rh allows for large differences in the ionization structure 
against small differences in rh and therefore that the '1 keV 
feature may not be ubiquitous in NLSyls. Indeed, a simi- 
lar feature is not present in NAB0205+024. A line with the 
same equivalent width as in PG1244+026 (60 eV) would 
have been detected in the NAB0205+024 SIS spectrum (the 
90 % upper limit is only 20 eV). In the accretion disc sce- 
nario this would imply a different ionization state (higher or 
lower) of the matter or a significant metal underabundance 
in NAB0205+024. 

5.3 Variability 

The mean 0.5-10 keV luminosity measured by ASCA in 
PG1244+026 and NAB0205+024 differs by an order of mag- 
nitude: 7.3 x 10 43 and 6.4 x 10 44 erg s _1 respectively (assum- 
ing isotropic emission, Ho = 50.0 and qo = 0.0). The opti- 
cal (3000 A) monochromatic luminosities differs even more: 
4 x 10 43 and 9.4 x 10 44 erg s _1 respectively. The variations 
seen in the NAB0205+024 light curve imply an efficiency in 
the conversion of matter into radiation greater than 1.6 % 
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(e.g. Fabian 1984). Taken at face value, this excludes ther- 
monuclear reactions as the origin of the observed X-ray lu- 
minosity for which the upper limit on the efficiency is in this 
case about 0.7 %. 

The similar variability observed in the two quasars 
agrees with the Fiore et al. (1998a) finding that the variabil- 
ity properties of (PG) quasars are correlated with the shape 
of the soft X-ray spectrum and the width of the Balmer lines, 
and so possibly then with the accretion rate, in the scheme 
of Pounds et al. (1995), and Laor et al. (1994), (1997). 



6 CONCLUSIONS 

ASCA observations of two steep soft X-ray quasars have 
shown that: 

(i) The X-ray continuum of the two quasars flattens by 
Act = 0.5 — 1 going toward high energies. Similar results 
were obtained by authors on some half dozen steep soft X- 
ray quasars. 

(ii) PG1244+026 shows a significant feature in the '1 keV 
region. Similar features were again reported in other steep 
soft X-ray quasars. The data are not good enough to dis- 
criminate between a broad emission line centered at 0.95 
keV (quasar frame) or an absorption edge at 1.17 keV, or 
an absorption notch at 1.22 keV. 

Line emission could be due to reflection from an highly 
ionized accretion disk, in line with the view that steep soft 
X-ray quasars are emitting close to the Eddington luminos- 
ity. Photoelectric edge absorption or resonant line absorp- 
tion could be produced by gas outflowing at a large veloc- 
ity (0.3-0.6 c). In these absorption models significant cold 
(i.e. oxygen less ionized than OVI) absorption in excess of 
the Galactic is required. This would imply an increase by 
a factor 2-3 of the cold column with respect to a previous 
PSPC observation or a peculiar ionization structure. In nei- 
ther the emission or absorption cases the SIS resolution is 
good enough to identify unambiguously the ions responsible 
for the feature. The high resolution and high throughput 
of the low energy gratings and spectrometers of AXAF and 
XMM are clearly needed to shed light on this puzzling case. 

(iii) The two quasars show similar variability properties 
(flux variations up to a factor of 2 in 10 ks) despite a fac- 
tor of ten difference in the X-ray observed luminosity. This 
agrees with the Fiore et al (1998a) finding that the variabil- 
ity properties of radio-quiet quasars are correlated with the 
shape of the X-ray spectrum, the width of the Balmer lines 
and so possibly with the accretion rate. 

F.F acknowledges support from NASA grants NAG 5- 
2476 and NAG 5-3039, B.J.W. acknowledges support from 
ASC contract NAS8-39073. 
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